AN EULERIAN SPACE-TIME FINITE ELEMENT METHOD FOR 
DIFFUSION PROBLEMS ON EVOLVING SURFACES 
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Abstract. In this paper, we study numerical methods for the solution of partial differential 

equations on evolving surfaces. The evolving hypersurface in W 1 defines a d-dimensional space- 

C^") , time manifold in the space-time continuum R ti + 1 . We derive and analyze a variational formulation 

for a class of diffusion problems on the space-time manifold. For this variational formulation new 

well-posedness and stability results are derived. The analysis is based on an inf-sup condition and 

£*\J ' involves some natural, but non-standard, (anisotropic) function spaces. Based on this formulation a 

, discrete in time variational formulation is introduced that is very suitable as a starting point for a 

^v discontinuous Galerkin (DG) space-time finite element discretization. This DG space-time method 

_4 , is explained and results of numerical experiments are presented that illustrate its properties. 

1. Introduction. Partial differential equations (PDEs) posed on evolving sur- 
faces arise in many applications. In fluid dynamics, the concentration of surface active 
agents attached to an interface between two phases of immiscible fluids is governed by 
a transport-diffusion equation on the interface |12| . Another example is the diffusion 
of trans-membrane receptors in the membrane of a deforming and moving cell, which 
is typically modeled by a parabolic PDE posed on an evolving surface [2]. 

Recently, several approaches for solving PDEs on evolving surfaces numerically 
have been introduced. The finite element method of Dziuk and Elliott [6] is based 
on the Lagrangian description of a surface evolution and benefits from a special in- 
variancc property of test functions along material trajectories. If one considers the 
Eulerian description of a surface evolution, e.g., based on the level set method [TT] . 
then the surface is usually defined implicitly. In this case, regular surface triangula- 
\f} • tions and material trajectories of points on the surface are not easily available. Hence, 

lO \ Eulerian numerical techniques for the discretization of PDEs on surfaces have been 

studied in the literature. In [TJ [T5] numerical approaches were introduced that are 
based on extensions of PDEs off a two-dimensional surface to a three-dimensional 
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neighbourhood of the surface. Then one can apply a standard finite element or (as 
was done in [Bill]) finite difference disretization to treat the extended equation in M 3 . 
The extension, however, leads to degenerate parabolic PDEs and requires the solution 
of equations in a higher dimensional domain. For a detailed discussion of this exten- 
sion approach we refer to [TTJ [3 13] . A related approach was developed in [8] , where 
S^ ' advection-diffusion equations are numerically solved on evolving diffuse interfaces. 

A different Eulerian technique for the numerical solution of an elliptic PDE posed 
on a hypersurface in M 3 was introduced in [T^ Q3] . The main idea of this method is to 
use finite element spaces that are induced by the volume triangulations (tetrahcdral 
decompositions) of a bulk domain in order to discretize a partial differential equation 
on the embedded surface. This method does not use an extension of the surface partial 
differential equation. It is instead based on a restriction (trace) of the outer finite 
element spaces to the (approximated) surface. This leads to discrete problems for 
which the number of degrees of freedom corresponds to the two-dimensional nature 
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of the surface problem, similar to the Lagrangian approach. At the same time, the 
method is essentially Eulerian as a surface is not tracked by a surface mesh and may 
be defined implicitly as the zero level of a level set function. For the discretization 
of the PDE on the surface this zero level then has to be reconstructed. Optimal 
discretization error bounds were proved in [15] . The approach was further developed 
in [4j [16], where adaptive and streamline diffusion variants of this surface finite element 
were introduced and analysed. These papers [HI [14] [4] [16] , however, treated elliptic 
and parabolic equations on stationary surfaces. 

The goal of this paper is to extend the approach from |15] to parabolic equa- 
tions on evolving surfaces. An evolving surface defines a three-dimensional space-time 
manifold in the space-time continuum R 4 . The surface finite element method that we 
introduce is based on the traces of outer space-time finite element functions on this 
manifold. The finite element functions are piecewise polynomials with respect to a 
volume mesh, consisting of four-dimensional prisms (4D prism = 3D tetrahedron x 
time interval). For this discretization technique, it is natural to start with a variational 
formulation of the differential problem on the space-time manifold. To our knowledge 
such a formulation has not been studied in the literature, yet. One new result of this 
paper is the derivation and analysis of a variational formulation for a class of diffusion 
problems on the space-time manifold. For this formulation we prove well-posedness 
and stability results. The analysis is based on an inf-sup condition and involves some 
natural, but non-standard, (anisotropic) function spaces. A second important result 
is the formulation and analysis of a discrete in time variational formulation that is 
very suitable as a starting point for a discontinuous Galerkin space-time finite ele- 
ment discretization. Further, we present a finite element method, which then results 
in discretization (in space and time) of a parabolic equation on an evolving surface. 

The discretization approach based on traces of an outer space-time finite element 
space studied here is also investigated in the recent preprint [TU]. In [TU] there is 
no analysis of the corresponding continuous variational formulation, which is the 
main topic of this paper. On the other hand, in |10| one finds more information on 
implementation aspects and an extensive numerical study of properties (accuracy and 
stability) of this method and some of its variants. We only very briefly comment 
on implementation aspects and illustrate accuracy and stability properties of the 
discretization method by results of a few numerical experiments. 

In this paper, we do not study discretization error bounds for the presented 
Eulerian space-time finite clement method. We expect that such error bounds can be 
derived using the tools developed for the analysis of the stationary case in [TS] HH US] 
and using the variational setting for the evolving surface case presented in this paper. 
This is a topic of current research. 

The remainder of the paper is organized as follows. In section[2]wc review surface 
transport-diffusion equations and introduce a space-time weak formulation. Some 
required results for surface functional spaces are proved in section [3] A general space- 
time variational formulation and corresponding wcll-poscdncss results are presented 
in section [4] A semi-discrete in time method is analyzed in section [5] A fully discrete 
space-time finite element method is considered in section [5] Section [7] contains results 
of some numerical experiments. 

2. Equation for diffusion on an evolving surface. Consider a surface T(t) 
passively advected by a smooth velocity field w = w(x, £), i.e. the normal velocity 
of r(i) is given by w ■ n, with n the unit normal on T(t). We assume that for all 
t G [0,T], r(i) is a smooth hypersurface that is closed (dT = 0), connected, oriented, 
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and contained in a fixed domain Q C K 3 . The conservation of a scalar quantity u 
with a diffusive flux on T(t) leads to the surface PDE (cf. [T5]): 

u + (div r w)u - v d A r u = on T(<), t € (0,T], (2.1) 

with initial condition u(x, 0) = Uq(x) for cc S T(0). Here u = §f + w • Vu denotes the 
advective material derivative, divr := tr ((/ — nn T )V) is the surface divergence and 
Ar is the Laplace-Beltrami operator, Vd > is the constant diffusion coefficient. 
Consider the space-time manifold 

S = |J T(t) x {t}, S C M 4 , 
te(o,T) 

and let H 1 (S) be the usual Sobolev space on S. The following weak formulation of 
(|2.1[) was shown to be well-posed in [6j: Find w € H 1 (S) such that w(-,0) = «o and 
for almost all t £ (0, T] 

/ to + uv div r w + y d V r u ■ V r v ds = for all v(-,t) e iJ 1 (r(t)). (2.2) 

Jr(t) 

Here Vr is the tangential gradient for T(t). The formulation (|2.2j) is a natural starting 
point for finite clement methods based on the Lagrangian description of the surface 
evolution. It is, however, less suitable for the Eulerian finite element method that 
we introduce in this paper. Our discretization method uses the framework of space- 
time finite element methods. Therefore, it is natural to consider a space-time weak 
formulation of (|2.ip as given below. We introduce the space 

H = {veL 2 (S) | ||VH| L 2 (5) <co} 

endowed with the scalar product 

(u,v) H = (U,V) L 2( S ) + (VrW,V r u)L2(5), (2-3) 

and consider the material derivative as a linear functional on H. The subspace of 
all functions v from H such that v defines a bounded linear functional form the trial 
space W. A precise definition of the space W is given in section 13.21 We consider 
the following weak formulation of (|2.1|) : Find u € W such that 

(u, v) + I / uv divrw + v^Vru ■ Vr« dsdt = for all v € H, 

Jo Jr(t) (2.4) 

u(-,0) = u . 

We shall derive certain density properties for the spaces W and H, which we need 
for proving the well-posedncss of (|2.4|) . Actually, we show well-poscdness of a slightly 
more general formulation, which includes a possibly nonzero source term and, instead 
of (divrw)w, a generic zero order term. Our finite element method will be based on 
(HID rather than ([2~2")h 

3. Preliminaries. In this section, we define the trial space W and prove that 
both the test space H and the trial space W are Hilbert spaces, and that smooth 
functions are dense in H and W. We also prove that a function from W has a well- 
defined trace as an element of L 2 (T(t)) for all [0, T]. In the setting of a space-time 
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manifold, the spaces H and W are natural ones. In the literature, however, we did 
not find any analysis of their properties. The necessary results are established with 
the help of a homeomorphisms between if, W and the following standard Bochner 
spaces H and W: 

H;=L 2 (0,T;H l (T )), W := { u G H \ -^ G H' }, T := T(0). (3.1) 

In the next subsection, we collect a few properties of the Bochner spaces H and W 
that we need in our analysis. 

3.1. Properties of the spaces H and W. The spaces H and W are endowed 
with the norms 

~ * „ , „ (, ,> .*.> ' 

and \\u\\ W =l\\u\\ s + \\ — \\ s , 

We start with the following well-known result. 

Lemma 3.1. The space C^°(r x (0,T)) is dense in H. 

Proof. The inclusion Cfi°(T x (0, T)) C H is trivial. By construction of the 
Bochner space, the set of simple functions { Yl7=i XB^i \ 4>i G ^(To), n G N} is 
dense in H; here {B>i} is any set of n mutually disjoint measurable subsets of (0,T). 

For XBii'i there exists gi G C^°(0,T) and (pi G C°°(ro) such that \\xBi4>i ~ Qi&Wji 1S 
arbitrary small. This completes the proof. D 

For w G H we define the weak time derivative through the functional 

(%<!>)-■=- j j T w{ ti% dsdt for ^eCgPWO.T)). (3.2) 

Then ^ G H' iff there is a constant c such that 




dw 

eF'« 



< 



c||<% for all ^^(r o x(0,T)). 



Remark 3.1. The definition of W in (|3.1|l . based on the weak time derivative 
is equivalent to the following more standard one: w G H is an element of W iff 
there exists z G L 2 (0,T; F _1 (r )) such that 

(z(i),v) H -i xH i<t>(t)dt = - / w(t)v4>' (t) ds dt (3.3) 

Jo Jt 

for all v G ^(To), G Cg°(0,T). The definition of V? in dSTTJ) is more convenient 
for the analysis that follows. 

We recall a few results for the space W. 

Lemma 3.2. The set 

e = \ f>V* | fer(r„), hgn! 

is dense in W . For u G W the function t — > u(t) = u(-,t) is continuous from [0,T] 
into L 2 (Tq). There is a constant c such that 

maxJu(t)|| L 2 (rn) < c\\u\\^ for all u G W. (3.4) 



Proof. Proofs are given in standard textbooks, e.g., [TI5] Proposition 23.23. The 
density result is usually proved with C°°(ro) replaced by -H^fTo) in the definition of 
C. The result with C°°(r ) follows from the density of C°°(r ) in H l (T ). U 

3.2. The spaces H and W. We assume that the space-time surface S is suf- 
ficiently smooth, cf. section 13.31 for a more precise description of the smoothness 
assumptions. Due to the identity 

/ f(a,t)dsdt= //(.s)(l + (wn) 2 )-5d S , (3.5) 

'o Jr(t) Js 

the scalar product L J r(t \ vw ds dt induces a norm that is equivalent to the standard 
norm on L 2 (S). Therefore, one can equivalently define the norm on H by 

r T 

2 ._ / / „,2 , |Y7_„,|2, 



M\h-= / / v 2 + \V r v\ 2 dsdt. (3.6) 

Jo Jr(t) 

The space H is a Hilbert space, and H =-> L 2 (S) =-> H' forms a Gelfand triple (cf. 
Lemma 13.51 below) . 

Recall the Leibniz formula 

/ v + vdiv T wds = — / vds, veC^S), (3.7) 

which implies the integration by parts identity: 

/ / iiv + vu + uv divpw ds dt 

h Jr ^ (3.8) 

= / u(s,T)v(s,T)ds- u(s,0)v(s,0)ds for all u,v € C^S). 

Based on (|3.8p we define the material derivative for u <G H as the functional it: 

(u,<f>} = - / M0 + u0div r wrfsrft for all <j) e C^{S). (3.9) 

Jo Jr(t) 

Assume that for some u € H the norm 

,, . ,, (u,d>) 

\\u\\h> = sup 



is bounded. In Lemma 13.51 we prove that Cq(S) is dense in H and therefore ii can 
be extended to a bounded linear functional on H. In this case, we write u G H' and 
define the space 

W = {ueH | iieH'}, with \\u\\ 2 w :=\\u\\ 2 H + \\u\\ 2 H ,. 

In section [BT^l we prove that W is a Hilbert space and C 1 (5) is dense in W. Note that 
the space W is larger than the standard Sobolev space iJ 1 (5). 

Remark 3.2. From the definition of the weak material derivative in Q3.9p and 
the density result of Lemma 13.51 it follows that for u £ C 1 (S) we have 

(u, v) = / uv ds dt for all v £ H. 

Jo Jv(t) 



3.3. Homeomorphism between {H, W} and {H, W}. To define a homeo- 
morphism we use a Langrangian mapping from the space-time cylinder T x [0, T] to 
the space-time manifold S, see also [5]. We assume that the velocity field w and T 
are sufficiently smooth such that for all y eTq the ODE system 

*(y,0)=V, -^{y,t)=w{${y,t),t), t€[0,T], 

has a unique solution x := &(y,t) 6 T(t) (recall that T(t) is transported with the 
velocity field w). The corresponding inverse mapping is given by $ _1 (a;, t) := y G Tq, 
x £ T(i). The Lagrangian mapping $ induces a bijection 

F: r o x[0,T]^S, F(y,t) := (${y,t),t). 

We assume this bijection to be a C 2 -diffcomorphism between these manifolds. 
For a function u defined on S we define u = u o F on To X (0, T): 

«(»,«) = «($(y,t),t) = «(a; > t). 

Vice versa, for a function u defined on T x (0, T) we define u = u o _F _1 on 5: 

u(ar,t) = u($ _1 (a;,t),£) = u{y,t). (3.10) 

By construction we have 

u(x,t) = —(y,t). (3.11) 

Now we prove that the mapping u — > u defines a linear homeomorphism between H 
and H , and also between W and W . 

Lemma 3.3. The linear mapping u — » u /rom (|3.10[) defines a homeomorphism 
between H and H . 

Proof. For any fixed t G [0,T], we obtain u(t) := u(-,t) G H 1 ^) iff u(t) = 
u{-,t) G H 1 (r(i)). Let u(t) G H 1 (T Q ) for all £ G [0,T]. Due to the smoothness 
assumptions on i* 1 , there are constants c\,cq > 0, independent of u and £, such that 

c ||u(t)|| ffl( r ) < \\u{t)\\ H i {m) < Cl ||u(t)|| ff i (ro ) for all t G [0,T]. (3.12) 

Hence, a £ H iff w € If holds, and the linear mapping u — > it is a homeomorphism 
between 7J and H . D 

For the further analysis, we need a surface integral transformation formula. For this 
we consider a local parametrization of Fo, denoted by \i : M. 2 — > Tq, which is at least 
C 1 smooth. Then, $ o /i := $(/x(-),i) defines a C 1 smooth parametrization of T(t). 
For the surface measures d J and ds on To and T(t), respectively, we have the relation 

^JVr*(.t)/i.xVr*(-,t) / g 

\HX X fj,y\ 

with [i x = /Ua;(jU _1 (-)), and similarly for /j y . Recall that Vr/ denotes the r(i)-surface 
gradient of a scalar function / defined on T(t) for any fixed t. Using this integral 
transformation formula, for u G H and G Cq(<S) we obtain 



/ / u<j> + u<f> divpW ds dt = — / / (u— — hM</>divrw)7(is"<i£ 
Jo Jr(t) Jo Jr dt 

(3.14) 



Lemma 3.4. The linear mapping u — > u from (|3.10[) defines a homeomorphism 
between W and W . 

Proof. The proof makes use of the formula (|3.14[) . Take u £ H with u £ W, and 
4> £ Cq{S). We use the notation ~ and < if the constants are independent of u and 
(f>. Due to the C 2 -smoothness assumption on F (and thus <f>) the function 7 defined 
in (|3~T3l) is ^-smooth on T x [0,T]. Define $ := 7 £ C^(T x (0,T)). Due to 
Lemma |3~31 we have ||V>|||j < ||<^||# — ll^llff- Therefore, we can estimate 



/ (u— — h M^divrw)7<is"rf<| 



-To 

r r Pi 1 r r f) r r 

^ / / u-^— d'sdtl + \ I / ucf)—-dsdt\ + \ I / u<f> divrw 7 ds dt\ 
Jo Jr a dt ' 'Jo Jt dt Jo Jr 

< \\u\\w\$Wh + MhUWh < MwMIh- 

Hence, u £ W and ||w||w ^5 ||u||tt? holds. With similar arguments one can show that 

if u £ W, then u £ W and ||m||^ < \W\\w holds. For this, instead of the surface 
transformation formula (|3.13[) one starts with the formula 

d's = 7= — =rz — n — ; — \ — - — J! , /, — w — \ — ; — r ds —: y(-,t) ds, (3.15) 

iVr^^-^-.O^Kx V r $($- 1 ('^)^KI 

with /i x = jU x (($ o /x) _1 (-)), and similarly for pL y . For u £ H and G Co(r x (0,T)) 
we have 

— ,4>) = - / / u—dsdt=- / / ucjrjdsdt 
dt I Jo Jr dt J J T(t) 

Now we note that 7 is C 1 -smooth on S. To check this, it is sufficient to show that 
the denominator in (|3.15p is uniformly bounded away from zero on S. For x £ T(t) 
and with x = ($ o yi){z) one can rewrite the denominator as 

|V r $($ _1 (z>*),*)/4r(.z) x Vr$($ _1 (£,i),*W-z)l 

= |Vr*(M(*).*)M*) x Vr*(j«(«),t)A»»(«)| / g 16) 

= l^( $ °M)(^)x t-($o M )(z)|. 

From the fact that $ o /1 is a C^-smooth paramctrization of T(t) it follows that the 
quantity on the right-hand side of (|3.16j) is uniformly bounded away from zero. Hence, 
the function 7 is C^-smooth and we can use the same arguments as above to derive 
ll^llw ^S IMIw- This implies that u — > u is a homeomorphism between W and W. □ 

3.4. Properties of H and W . The homeomorphism established in i j3.3l helps 
us to derive density results for the spaces H and W and a trace property similar to 
the one in (|3.4[) . 

Lemma 3.5. H is a Hilbert space. The space Cq(S) is dense in H. The spaces 
{H,L 2 (S),H'} form a Gelfand triple. 

Proof. Let L : H — > H denote the mapping given in p.lO|) . Since L is a linear 
homeomorphism between, the space H is complete and so this is a Hilbert space. For 
<fi £ Cq(Fo x (0,T)) C H we have, due to the smoothness assumptions on F, that 
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4> = Ltfi G C 1 (5). Furthermore, from 4>(x,t) = 0($ 1 (x,t),t) it follows that </> has 
compact support. Hence, 4> G Cq(S). From this we get i[Co(r x (0,T))] C C^(S). 
Since Cq(Fo x (0,T)) is dense in H and L is a homcomorphism, this implies that 
Cq(S) is dense in H. Since Cq(5) is also dense in L 2 (S), the space H is dense in 
L 2 {S). Hence, {H,L 2 (S),H / } is a Gelfand triple. □ 

For £ G [0,T] and it G C l (S) denote by u — > U|p(t) a trace operator: U|r(t)(x) = 
it(a:,£), x € r(i). In section [5] we analyze a discontinuous Galerkin method in time. 
For such a method, one needs well defined traces of this type. For a smooth function 
u(x,t) defined on the cylinder To x [0, T], it is obvious that one can define at any 
time t G [0,T), the right limit «+(•) = lim u(-,t + 6) on To. Similarly a left limit 

(5— >+0 

function u- is defined for t G (0,T]. For a sufficiently smooth function u on S, due 
to the fact that the domain T(t) where the trace has to be defined varies with t, it is 
less straightforward to construct such left and right limit functions. To this end, for 
u G C X {S) and a given t G [0, T] we define u s : T(t) -> R by 

UiO,*) :=u(F($- 1 (-,i),t + ,5)), <5 such that i + <5 e [0,T]. (3.17) 

Note that w^ = uir(t) holds when 5 = 0. Right and left limits on T(t) are defined as 

u+(;t)= lim u s (;t), forte [0,T), 

(3.18) 
u_ (-,*)= lim w 5 (-,t) foriG(0,T]. 

(5— > — 

Below we show that for functions from W the trace and one-sided limits arc well- 
defined and can be considered as elements of L 2 (T(t)). 

The next theorem gives several important properties for our trial space. 

Theorem 3.6. W is a Hilbert space and has the following properties: 
(i) C 1 (iS) is dense in W . 

(ii) For every t G [0, T] the trace operator u — > W|r(t) can be extended to a bounded 
linear operator from W to L 2 (T(t)). Moreover, the inequality 



max 

0<t<T 



ll«|r(t)llL 2 (r(t)) < c\\u\\ w for all u G W, (3.19) 



holds with a constant c independent of u. 
(Hi) Take t G [0, T) and let e > be sufficiently small, such that t + e < T. For any 
u G W the mapping 5 — > Ug(-,t) defined in (|3.17|1 is continuous from [0,e] 
into L 2 (T(t)). The same assertion is true for t G (0,T] and suitable e < 0. 
For u G W the one-sided limits (|3.I8[) are well-defined. 
Proof. Since the mapping L : u — > u given by (|3.10[) is a linear homeomorphism 
between W and W, the space W is complete and so this is a Hilbert space. 

(i) Let C be the set as in Lemma 13.21 which is dense in W . One easily checks 
L{C) C C 1 (iS). Since L(C) is dense in W, this implies that C 1 (5) is dense in W. 

(ii) Take u G C 1 ^) and define U|r(t) := u(-,t) G L 2 (T(t)). Using (|3T3]) . 
Lemma 13^21 and Lemma l3T4l we get 



lhr(t)IU 2 (r(t)) < c||u(t)|| L 2 (ro) < c||u||^ < c||u|| w , 

where the constant c can be assumed to be independent of t due to the smoothness 
of 7 in (|3.I3[) . From this, the result in (|3.I9I) follows by a density argument. 



(iii) Take a fixed t G [0, T) and sufficiently small e > 0. Take <$i,<$2 S [0,e], For 
x G r(i) we use the substitution j/ := $ _1 (.t,£) G r and the integral transformation 
formula as in the proof of Lemma 13.31 resulting in: 



\\u s d-,t)-u S2 (-,t)\\ L 2 {r{t)) = \\u(F($-\-,t),t + 6 1 ))-u(F($- 1 (-,t),t + 6 2 ))\\ L 2 (m) 

< c\\u(-,t + St) -u(-,t + 6 2 )\\ L 2 ( r o) , 

with a constant c independent of t. Hence, the continuity of the mapping S — > us{-, t) 
follows from the continuity result for u in Lemma 13.21 Due to the continuity of the 
mappings, the one-sided limits are well-defined. D 

Corollary 3.7. For all u, v G W, the integration by parts identity holds: 

(u,v) + (v,u) + / / uv divrw ds dt 

Jo Jt(X) (320) 

u(s,T)v(s,T) ds — / u(s, 0)v(s, 0) ds. 

T(T) Jr(o) 

Proof. Follows from the identity (|3.8p and Theorem 13.61 □ 

4. Well-posedness of weak formulation. Using the properties of H and W 
derived above, we prove a well-posedness result for the weak space-time formulation 
(|2.4p of the surface transport-diffusion equation (|2.1|) . As usual, we first transform 
the problem (|2.1[) to ensure that the initial condition is homogeneous. To this end, 
consider the decomposition of the solution u = u + u°, where u° : S — >• K is chosen 
sufficiently smooth and such that 

u°(x, 0) = Uq(x) on To, and — / u° ds = 0. 

"* Jr(t) 

One can set, e.g., u° = (wo 0< l )_1 )(7 -F 1_1 )~ 1 , with 7 from p,13|) . Since the solution of 
(|2.1[) has the mass conservation property ^ J r , t -, u ds = 0, and § T i \ u° ds = f r , \ u ds 
by the choice of u°, the new unknown function u satisfies u(-, 0) = on Tq and 

/ uds = for alHG [0,T]. (4.1) 

Jr(t) 

For this transformed function the surface diffusion equation takes the form 

u + (divrw)u - v d Aru = f on T(t), t G (0, T], 
u(-,0) = on r . 

The right-hand side is now non-zero: / := — u° — (divrw)u + i^Ar - " . Using the 
Leibniz formula (|3.7[) and the integration by parts over T(t), one immediately finds 
f m fds = foralltG [0,T]. 

In the analysis below, instead of the (transformed) surface diffusion problem (|4.2p 
we consider the following slightly more general surface PDE: 

u + au- Vd&ru = f on T(i), t G (0, T], 
u = on To, 

with a G L°°(S) and a generic right-hand side / G L 2 (S), not necessarily satisfying 
the zero integral condition. We use the notation a^ := \\a\\i,<x^sy 
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Wc define the inner product and symmetric bilinear form 

(u,v)o= / / uvdsdt, a(u, v) — ^(Vru, Vrf)o + { au : v)o, u,v £ H. 
Jo Jr(t) 

This bilinear form is continuous on H x H: 

a(u,v) < (v<i + oc^WuWhWvWh- (4.4) 

Consider the subspace of W of all function vanishing for t = 0: 

W:={v£ W | w(-,0)=0 onTo}. 

o 

The space W is well-defined, since functions from W have well-defined traces on T(t) 
for any £ £ [0,T], see Theorem 13.61 The weak space-time formulation of (|4.3p reads: 

o 

Given f & H', find u £W such that 

(u, u) + a(u, v) = (/, v) for all t> e H. (4.5) 

In the remainder of this section we prove that this variational problem is well- 
posed. The analysis is based on the continuity and inf-sup conditions, cf. [9]. The 
continuity property is straightforward: 

\ (ii,v} + a(u,v)\ < (1 + v,i + c(oo)\\u\\w\\v\\h for all u £ W, v £ H . 



The next two lemmas are crucial for proving the well-posedness of (|4.5|) . 
Lemma 4.1. The inf-sup inequality 

■ t (u,v) +a(u,v) 

mf sup > c s (4.6) 

holds with some c s > 0. 

o o 

Proof. Take u £ W . In (|4.5[) we take a test function v = u 1 := e rt u £ W, with 
7 := 2{vd + ||a — g divrw|j^=o(5)). We note the identity: 

{u-f,u) = (it, u-y) - j(u, u 7 ) . (4.7) 

From (|4~71) . (|X20|) . and condition u(0) = 0, wc infer 

(w,u 7 ) = -((u,u y ) + (v.j,u)) + -7(w,%)o > --(-u,it 7 div r w) + -7(m,u 7 ) . 

This and the choice of 7 implies (u, tt 7 ) + (qm, m 7 )o > ^d(u,u 7 )o > ^de~ 7T ||ii||Q. 
Combining this with (Vrw, Vru 7 )o > e ~ 7T ||^rM||o, we S e t 

(it, Uj) + a(u, Uj) > v^e ' \\u\\ H . (4-8) 

This establishes the control of \\u\\h on the right-hand side of the inf-sup inequality. 
We also need control of ||u||ff' to bound the full norm ||u||w- This is achieved by 
using a duality argument between the Hilbcrt spaces H and H'. 

By Riesz' representation theorem, there is a unique z £ H such that (ii, v) = 
(z,v)h for all v £ H, and \\z\\h = \\u\\h' holds. Thus wc obtain 

(u,z) = (z,z)h = \\u\\h>- 
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Therefore, with the help of (|4.4[) . we get 

■■)>\\4 2 H-^\\U\\ 2 H -I\\z\\l 

(4.9) 



(u, z) + a(u, z) = \\z\\% + a(u, z) > \\z\\ 2 H - -^||u||h- ^\\z\\ 2 h 



= 2 NIh' ~ yll^llff' with C l = ( V d + "oo) 2 . 

This establishes control of ||u||h' at the expense of the ff-norm, which is controlled 
in (|4.8|) . Therefore, we make the ansatz v = z + ^ u 7 G H for some sufficiently large 
parameter /j, > 1. We have the estimate 

IMI* < IMIh + mIKIIh < II«IIh' + mIMIh < /"V2|klk- (4.10) 

From dUSJl, dH}, and (|^7T0|) we conclude 

(ii,v) +a(u,v) > -\\u\\h' + {vvde~ lT - -^)\\u\\ 2 H . 
Taking \x := ^ (a + l)e 7T , we get 

1 /2 

(u,w) +a(u,v) > -\\u\\ 2 w > —h^WuWwWvWh- 

This completes the proof. D 

Remark 4.1. A closer look at the proof reveals that the stability constant c s in 
the inf-sup condition (|4.6[) can be taken as 

c s = ^±(l + Vd + aoo )- 2 e- 2T ^+~ c \ c=\\a-\ div r w|| L ^ (5) . 

This stability constant deteriorates if V& \. or T —*■ oo. We do not consider the 
singularly perturbed case of vanishing diffusion. Without assumptions on the sign of 
a and the size of the velocity field w (which is part of the material derivative) , there 
may be (exponentially) growing components in the solution and thus the exponential 
decrease of the stability constant as a function of T can not be avoided. In special 
cases, the behavior of the stability constant may be better, e.g. bounded away from 
zero uniformly in T. We comment on this further after Theorem 14.31 

o 

Lemma 4.2. If (it, v) + a(u, v) = for some v £ U and all u £ W , then v = 0. 

o 

Proof. Take v £ H such that (u, v) = —a(u, v) for all u &W. For all u £ Cq(S) C 

o 

W we have by definition 

(ii,u) = — / viidsdt — / / uv divrw ds dt 

Jo Jr(t) Jo Jr(t) 

= —(u,v)— / / uvdivr'wdsdt = a(u, v) — (u, vdivrw)o. 
Jo Jr(t) 

Since the functional u — > a(u,v) — (u, v divrw)o is in H' , we conclude v £ H' , and 
thus v € W holds. From a density argument it follows that 

(v, u) = a(u, v) — (u, i>divrw)o for all u £ H (4-11) 

n 



holds. For all u G W C H we get 

(v, u) = a(u, v) — (u, v divrw)o = — (u, v) — (u, v divrw)o- 
This and (f3T20|) yield 



= (v,u) + (u.v) + (M,wdivr"w)o = / u(s,T)v(s,T) ds for all u £ W. 

Jr(T) 

This implies that v(-,T) = on T(T). We proceed as in the first step of the proof of 
Lemma RTTl We take in (|4.11[) u = v 1 = e^' yt v, with 7 := — 2(l+||a— i divrw|| L oc( 5 ) < 
0, and use (|4.7p . We obtain 



((v,v-y) + (v-y,v)) + -7(u,u 7 )o-a(u 7 ,u) + (v 7 ,vdiv r w) 



= (i), v-y) — a(u 7 , v) + (v 7 ,v divrw)o 

= i^{{i^i) + (' b 11 v)) + 2 
1 , ,1 



j{v,v 1 )o + -(u 7 ,udiv r w)o-a(u 7 ,u) 
< -0, v 7 ) - v d {\J T v- n V r w) < -(Ho + »d\\ Vrullo)- 
We conclude u = 0. D 

As a direct consequence of the preceding two lemmas we obtain the following 
wcll-posedness result. 

o 

Theorem 4.3. For any / £ H' , the problem (|4.5[) has a unique solution u £ W . 
This solution satisfies the a-priori estimate 

\\u\\w<c; l \\f\\ H '. 

Wc consider two special cases in which the inf-sup stability constant c s can be shown 
to be bounded away from zero uniformly in T, cf. Remark 14.11 
Proposition 4.4. Assume that there is c > such that 

a — — divrw > cq on S (4-12) 

holds. Then the inf-sup property (14.611 holds with c, = """t^coi _ 

Proof. Wc follow the arguments as in the proof of Lemma 14.11 Instead of v = u 7 
we take v = u as a test function. This yields 

(u,u) + a(u,u) > c\\u\\ H , c :— min-jV^, Cq}. (4-13) 

We set v = z + [iu [z as in the proof of Lemma l4~Tj) . Taking /1 := ^ti > 1 we get 

/• \ t \ ^ ll ll 2 v2 _ln 11 11 11 1^11^1^: Co} n n || n 

{u,v)+a{u,v) > -\\u\\ w > —fi \\u\\w\Hh = ,-,_, , ; rdMMMIff- 

2 4 V2(l + Vd + a oo y 

This completes the proof. □ 

As a second special case, we consider the surface diffusion problem (|4.2|) . A 
smooth solution to this problem satisfies the zero average condition (|4.1[) . Functions 
u from iJ satisfying (|4.1[) obey the Fricdrichs inequality 



/ \\7 r u\ 2 ds>c F {t) I u 2 ds for allie [0,T], (4.14) 
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with cp(t) > 0. The Fricdrichs inequality helps us to get additional control on the L 2 
-norm of u in proving the inf-sup inequality and so to improve the stability constant 
c s . Below we shall make more precise when a solution to (14.51) satisfies a zero average 
condition, and how the stability estimate is improved. 

o 

We introduce the following subspace of W: 

W:={ueW\ [ u(-,t)ds = for all i£ [0,T]}. 
Jr(t) 

Proposition 4.5. Assume a = divrw, / is sufficiently regular, e.g. f £ L 2 (S), 

o 

and Jw t \ f ds = for almost all t £ [0,T]. Then the solution u £ W of (|4.5|l belongs 
to W . Additionally assume that there exists a cq > such that 

div r w(x,t) + v d c F {t) > c Q for all x £ T(t), t £ [0,T] (4.15) 

holds. Then the inf-sup property (|4.6|l holds, with c s = """t^^oi an g yj re pi acec [ 
by the subspace W . 

o 

Proof. Let u £ W be the solution of (|4.5|) . Dchnc U(t) := J" r(t , u(-,i) <is. Using 
Theorem 13.61 (ii) we get 

U(t) 2 dt< f ||tt|r( t )||i a( r(t))|r(t)|«ft<c||«||^. 

o Jo 

Hence, U £ L 2 (0,T) holds. Take cf> £ C^(0,T), and thus £ Cg(«S). Note that 

U(t)<j)'(t)dt= - / ucf)' (t) dsdt = - \ I ucpdsdt 

o Jo Jr(*) Jo Jr(t) 

= (w, 0) + (u, 0div r w) o < c||«||w-||^||i? < c||«||w-||^||l2(o,t)- 

This implies that U £ H 1 ^, T) holds. Using J r{t) f(-,t) ds = for almost all < £ [0, T], 
we obtain, for arbitrary cfi £ Cq(0,T): 

T pT 

U'(j)dt = — I Ucf) 1 dt = (ii, cp) + (u, (^divrw)o = (u, </>) + a(u, cf>) 
o Jo 

= (/,0)o= / ^ / fdsdt = 0. 
jo Jr(t) 

Thus, [/ is a constant function. From f/(0) = J r u(-, 0) ds ~ it follows that U(t) = 

for all t £ [0, T] and thus u £ Wo holds. 

For the proof of the^inf-sup property, we follow the arguments as in the proof of 
Lemma l4.il Take u £ Wq. Instead of v = u 7 we take v = iiasa test function. Using 
the Fricdrichs inequality (|4.14[) and the assumption in (|4.15[) . we get 

(u,u) +a(u,u) > -(div r w,w 2 ) +^d||V r w||o > -(div r w + ^ d c F , w 2 ) + -^H^r^Ho 

> 2 C oll w llo+ 9^ll V r u llo > c\\u\\ 2 H , c:= - min{i/ d ,c }. 
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Take v :— z + jiu (z as in the proof of Lemma r4.il) . Taking jjl := ^ti > 1 we get 

(u,v) + a(u,v) > -\\u\\w > -j-it~ 1 \\u\\w\\v\\H = ttthTTTi ' ,° — ^TlMklMk- 

2 4 2V2(1 + Od + ctoor) 

This completes the proof. □ 

5. Time-discontinuous weak formulation. In this section, we study a time- 
discontinuous variant of the weak formulation in (|4.5[) . This new variational formula- 
tion is even weaker than (|4.5j) . However, it can be seen as a time-stepping discretiza- 
tion method for (|4.5j) and is better suited for the discontinuous Galcrkin discretization 
framework. 

Consider a partitioning of the time interval: = to < t\ < . .. < tjy = T, with 
a uniform time step At = T/N. The assumption of a uniform time step is made to 
simplify the presentation, but is not essential for the results derived below. A time 
interval is denoted by /„ := (t n _i, t„]. The symbol S n denotes the space-time interface 
corresponding to I n , i.e., S n := U t £i n T(t), and S := Ui<„<jv<5™- We introduce the 
following subspaccs of H: 

H n :={«£ H | v = on 5\5"}. 

For u £ H we use the notation u„ := u\ S n G H n . Corresponding to the space H n , we 
define a material weak derivative as in section f3. 21 For u € H n : 



(u,4>)j = — / / utj) + u(j) div r w ds dt for all 4>eCl{S n ). 
Jt n -i Jr(t) 

If for u £ H n the norm 



\\U\\H' n = SUp —-— 

4>eCi(S«) \\<P\\h 

is bounded, then by a density argument, cf. Lemma 13.51 u can be extended to a 
bounded linear functional on H n . We define the spaces 

W n = { V € H n | V € H' n }, \\ V f Wn = \\vf H + 11*11^ . 

Finally, we define the broken space 

N 

W b := ®% =1 W n , with norm \\v\\ 2 wb = £ \\v\\^ n . 

n=l 

Note the embeddings W C W b C H , and furthermore: 

(ii,v) = (ii n ,v) In for ueW,v E H n . (5.1) 

For u <E W n , we define the one-sided limits u™ = u + (-,t n ) and u™ = M_(-,t„) based 
on (|3.18[) . The limits are well-defined in L 2 (r(t„)) thanks to Theorem 13. 61 (item (hi)). 
At to and tjy only w+ and u 1 ^ are defined. For v € W b , a jump operator is defined by 
[v] n = vl - v™ G L 2 (r(t n )), n=l,...,N-l. For n = 0, we define [v}° = < . 

On the cross sections T(t„), < n < N, of S the L 2 scalar product is denoted by 

C0>0)t„ : = / ip(f>ds. 

JT{t n ) 
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In addition to a(-, •), we define on the broken space W b the following bilinear forms: 

N N 

d(u,v) = ^2d n (u,v), d n (u,v) = ([u] n -' L ,v r l- 1 ) tn _ 1 , (u,v) b = ^2{u n ,v n ) In 

n— 1 n— 1 

Instead of (|4.5[) . we now consider the following variational problem in the broken 
space: Given / € H' C {W b )' ', find u £ W b such that 

(u,v) b + a{u,v)+d(u,v) = (f,v) for all v e IF 6 . (5.2) 

We note that we use W b for the test space, since the term d(u, v) is not well-defined 
for an arbitrary v £ H . From an algorithmic point of view this formulation has the 
advantage that due to the use of the broken space W b = $* =1 W„ it can be solved in a 
time stepping manner. The final discrete method (section[6|) is obtained by combining 
the variational formulation (|5.2[) and a Galerkin approach in which the space W b is 
replaced by a finite clement subspacc. Before we turn to such a discretization method 
we first study consistency and stability of the weak formulation (|5.2[) . 

o 

Lemma 5.1. Let u £ W be the solution of (|4.5|) . Then u solves the variational 
problem (|5.2p . 



Proof. Let u £ W be the solution of (|4.5|) . Take a test function v £ W n . Since 

o 

for u £ W and for any t £ [0, T] the mapping (|3.17l) is continuous as a mapping from 
a sufficiently small interval [—£,0] or [0, e] to L 2 (T(t)), it follows that [u] n = for all 
n = 0, . . . , N — 1. Hence we get d(w, u) = 0. From (14.5|) and (|5.1[) we get 

(u n ,v)j + a(u n ,v) = (/, i>) for all w G IF™. 

and thus («, w) b + a(u, v) + d(u, v) = (/, v) holds. D 

For the stability analysis, we introduce the following norm on W b : 

1 1 N 

INI 2 := MAh + ^ III + 2 E IIM"" 1 !!?™-!- w e ^ 

n=l 

Note that W b is not a Banach space with respect to this norm. Stability of the 
time-discontinuous formulation is based on the next lemma. 

Lemma 5.2. Set 7 := 2(z^ + \\a — \ divrw||/,oo(s)). The following inf-sup in- 
equality holds: 

inf sud (^ v )b + d(u,v) + a(u,v) r J 

1111 sup ^_ e- i/j . 

¥«e«" o^eff' III"! NIj? 

Proof. We follow the arguments as in the first part of the proof of Lemma 14.11 
Let u = Yln=i Un ^ W b be given. As a test function we use v = u 7 = e ''u £ W b , 
and let w 7 . n = e~ Jt u n £ W n . From the definition of the weak material derivative we 
get 

and using (|3.20j) and the choice of 7, this yields 

AT 

1 

(U,U~,), = 



nib 

^ n=l 

JV 



9 X/ ( ("'i' M 7,™)/„ + ( u ~/,n,U n ) In +7(u„,U 7: „)oj 
n=l 

2 E ((«->«?)*« - «~ 1 ,<^ 1 )t„-i) - -(M,u 7 div r w) n + -7(u,u 7 ) . 
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Setting 7 := 2(z^ + \\a — \ divrw|| L oo(s)), wc obtain 

1 N 
(w,%) 6 + (a U ,u 7 )o>2E( e_7tn ll u -llL- e_7 '"" 1 IK" 1 |lL_ 1 )+^(^%)o 



n=l 

N 



> jE (e-^IKIIl -e-^-Mltir'llL-i) +W" 7r NS- 



(5.3) 



We also have (with w° := 0): 

AT W 

n— 1 n— 1 

i N 

= ^E^^diM-^iiL, + ktil, - iitOiu 

ra=l 

e 7 *"||u™ll 2 - ^-^'"-ill,,"- 1 !! 2 



(5.4) 



-5E( 



t„ ° 11"+ llt„-i 

AT 



+ V 7T ii^iil + ^E^ 7t "- 1 HN"- 1 ii^ 1 - 

Combining the results in (|5.3j) . ()5.4|) with (Vr"M, Vrw 7 ) > e _7T |)VrM||o we obtain 

(u, u 7 ) fo + d(w, u 7 ) + a(u, u 7 ) 

1 1 N 

> i/de-^llfiHi, + ie-^|| U w ||| + \ ^ e-^-IIN- 1 !!^-, > ^ T \\uf. 

n=l 

Finally note that ||ii 7 ||ij < \\u\\h < v d 2 \u\ holds. D 

We conclude that the weak formulations in (|4.5|l and (|5.2|l are equivalent in the 
sense that the unique solution of the former is the unique solution of the latter. Note 
that the norm controlled by the ellipticity condition in Theorem 15.21 is weaker than 
the one in Lemma 14.11 since the norm ||| ■ ||| does not contain the material derivative 
term ||w||_ff' that is part of the norm |H|w- As a direct consequence of lemmas 15.11 
and 15.21 wc obtain the following theorem. 

Theorem 5.3. For any f £ H' , the problem (|5.2|) has a unique solution u £ W b 
satisfying the a priori estimate 

H|]<e^||/||H<, 
with 7 defined in Lemma \5.2[ 

In the previous section, we noted that for the surface diffusion equation, which 
is a special special case of the surface equation (|4.3[) . the stability constant can be 
improved, cf. Proposition 14.51 An analogon holds for the discontinuous time-space 
problem (|5.2[) . 

Proposition 5.4. Assume a — divrw, / is sufficiently smooth, J r(t \ f ds = 
for almost all t £ [0, T], and there exists a cq > such that 

divrw(x,t) + VdCF(t) > cq for all x £ T(t), t £ [0,T]. 
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Then the unique solution u £ W b of the problem, (J5.2I) satisfies the a priori estimate 

|u|<2i/-»max{l > cb 1 }||/||^. (5.5) 

Proof. Thanks to assumptions and Proposition ^. 51 we know that the solution u to 
(|4.5[) is from W . Thus, u satisfies zero average condition (|4.1[) and due to Lemma [S~T1 
and Theorem 15.31 this is also the unique solution to (|5.2[) . Therefore, we may make 
use of the Fricdrichs inequality (|4.14p and prove the a priori estimate (|5.5[) following 
the lines of the proofs of Lemma T5. 21 with 7 = 0. D 

6. Finite element method. We introduce a conforming finite element method 
with respect to the time-discontinuous formulation (|5.2[) . The method extends the 
Eulcrian finite clement approach from |15[ 1141 I16j and uses traces of volume space- 
time finite element functions on S (the practical implementation uses a piccewisc 
linear approximation of S, as explained below). 

To define our finite element space Wh,At C W b , consider the partitioning of the 
space-time volume domain Q = fl x (0,T] c M. d+1 into time slabs Q n := il x I n . 
Corresponding to each time interval /„ := (t n -i,t n ] we assume a given shape regular 
simplicial triangulation T n of the spatial domain Q. The corresponding spatial mesh 
size parameter is denoted by h. For convenience we use a uniform time step At = 
tn — *n— l- Then Qh = 1J T n x In is a subdivision of Q into space-time prismatic 

n=l,...,N 

nonintersecting elements. We shall call Qh a space-time triangulation of Q. Note 
that this triangulation is not necessarily fitted to the surface S. Wc allow T n to 
vary with n (in practice, during time integration one may wish to adapt the space 
triangulation depending on the changing local geometric properties of the surface) 
and so the elements of Qh may not match at t = t n . 

For any n£ {1, . . . , N}, let V n be the finite element space of continuous piecewisc 
linear functions on T n . In section [51 wc comment on the case of higher order finite 
element spaces. First we define the volume time-space finite element space: 

Vh.&t •= {v : Q -*■ K I v(x,t) = 4>o(x) + t(j)\[x) on every Q n , with <p , (pi £ V n }. 

Thus, Vh,At is a space of piecewise bilinear functions with respect to Q„, continuous 
in space and discontinuous in time. Now wc define our surface finite clement space as 
a space of traces of functions from Vh,At on S: 

W h ,At ■= {w : S^R I w = v\s, veV h ,At}- 

The finite clement method reads: Given / £ H' . find Uh £ Wh,At such that 

(uh,Vh) b + a(uh,Vh)+d(uh,Vh) = (f,Vh) for all v h £ W h ,At- (6.1) 

As usual in time-DG methods, the initial condition for Uh(-,0) is treated in a weak 
sense and is a part of the right-hand side functional. Since, it/, £ C 1 (Q„) for all 
n = 1, . . . , N, the first term in (j6.ll) can be written as 



N rt n 

[Uh,V h ) b = v 



y~] / (~aT +w-'Vuh)v h dsdt. 



Note that this formulation allows one to solve the space-time problem in a time 
marching way, time slab by time slab. 

17 



At this moment, we have no proof of an inf-sup stability result for the fully 
discrete problem (|6.1[) . similar to the one in Lemma I5T2"] A specific technical difficulty 
in adopting the proof of Lcmma l5T2"l to the discrete case is that an exponentially scaled 
trial function u 7 = e 7 *w^ does not belong to the space of test functions anymore. One 
way out would be to modify the discrete space of test functions accordingly. Then, 
the required stability result easily follows, but this leads to a finite element method 
dependent on the parameter 7. This is not what we use in practice. 

Clearly, in the special case of condition (|4.12[) . one can prove the inf-sup inequality 
in the discrete case by following the arguments of Lemma 15.21 with 7 = and using 
the condition (|4.12[) to control the 7J-norm of a trial function, cf. Proposition 14.41 
Since for 7 = the test function is taken the same as the trial function, the inf-sup 
inequality becomes the ellipticity result 

(uh, u h ) b +a(u h , Uh) + d(v,h,Uh) > min{l, Co}|)ti/,|| 2 for all u h € Wh.At- (6.2) 



The space Wh.At has a finite dimension, and hence the ellipticity result (|6.2j) is suffi- 
cient for existence of a unique solution. We summarize this in the form of the following 
proposition. 

PROPOSITION 6.1. Assume (|4.12[) , then for any f <G H' the problem (|6.1[) has a 
unique solution Uh € Wh,At satisfying the a priori estimate 

|tt A |<i/-'max{l,Cb 1 }||/||^. 



The special case of the diffusion surface problem as described in propositions 14.51 
and 15.41 is less straightforward to handle, since in the discrete setting, the method is 
not pointwise conservative, i.e. the zero average condition (|4.1[) can be satisfied only 
approximately. Stability analysis of the finite clement problem (|6.1|) in this interesting 
case is the subject of current studies and results will be reported elsewhere. 

Before presenting numerical results, we comment on a few implementation aspects 
of our surface finite clement method. More details are found in the recent report [10] . 

By choosing the test functions Vh in ()6.1|) per time slab, as in standard space-time 
DG methods, one obtains an implicit time stepping algorithm. Two main implemen- 
tation issues are the approximation of the space-time integrals in the bilinear form 
(ii} ll V} l ) b + a(uh,Vh) and the representation of the the finite element trace functions 
in Wh At- To approximate the integrals, we first make use of the transformation for- 
mula (|3.5[) converting space-time integrals to surface integrals over S, and next we 
approximate S by a 'discrete' surface S . In our implementation, the approximate 
surface S h is the zero level of 4>h € Wt £ t , where <ph is the nodal interpolant of a 
level set function <fi(x,t), the zero level of which is the surface S. In the experiments 
considered in the next section, for <j> we use the signed distance function of S. To 
reduce the "geometric error" the interpolation <f>h £ W~ h £. can be done in a finite 

element space with mesh size h < h, At < At. In the examples considered in the 
next section, we use h = ^h, At = \At (one regular refinement of the given outer 
space-time mesh) 

For the representation of the finite element functions in Wh,At, we consider traces 
of the standard nodal basis functions in the volume space Vh,At- Obviously, only 
nodal functions corresponding to elements P £ Q such that P n S h 7^ should be 
taken into account. In general, these trace functions form a frame in Wh,At- A 
finite clement surface solution is represented as a linear combination of the elements 
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from this frame. Linear systems resulting in every time step may have more than 
one solution, but every solution yields the same trace function, which is the unique 
solution of (j6.ll) . 

7. Numerical experiment. In this section, we present results of a few nu- 
merical experiments to illustrate properties of the space-time finite element method 
introduced in Section [6) 

Example 1. First, we consider a shrinking sphere, represented as the zero level 
r(t) of the level set function <fi(x,t) = x\ + x\ + x 2 — 1.5 2 e~*. The initial sphere 
has a radius of 1.5 at t = 0. The corresponding velocity field is given by w = 
— |e~*/ 2 n, where n is the unit outward normal on T(t). Hence divpw = —1. Hence, 
the coefficient a is negative and the condition (|4.12[) is not satisfied. We choose a 
solution u(x,t) = (1 + xia;2a;3)e' and thus the right-hand side is given by f{x,t) = 
( — 1.5e* + ^-e 2t )xiX2Xz- The problem is solved by the space-time DG method in the 
time interval < t < 1. 

The outer domain is chosen as fi = [—2,2] x [—2,2]. The outer spatial triangu- 
lation T is a uniform tetrahedral triangulation of f2 with mesh size h = hk = 2~ k , 
k = 1, ...,4. This triangulation is chosen independent of n. The time step is taken 
as At = Ate = 2~~ e , (■ — 0, ... ,4. The outer space-time finite element space 14, a* 
and the induced trace space are defined as explained in the previous section. The 
smooth space-time manifold S = U te [o ! i]r(i) is approximated as follows. To a given 
outer space-time mesh one regular refinement (in space and time) is applied. The 
approximate, piecewise affine, surface S h = U^^iiT^t) is defined as the zero level 
set of the nodal interpolant of the level set function <fi(t) on this refined outer mesh. 
This approximate space-time manifold is constructed per time step. 

For the computation of discretization errors the continuous solution u is extended 
by constant values in normal direction. This extension is denoted by u e . The initial 
condition Uh(to) is the trace of the nodal interpolant to u e (to). We compute two 
approximate discrete errors as follows. 



err L 2 {tN) = \\u - u h \\ L 2 (Th{tN))) 



and 



At Ar_1 

err LHHl) = { — ||V r > e - « e )||| 2( r, l(to )) + Y. Ai H V r,K - »h)ll 



-l !/ 2 

Vr h (u e -u h )\\ L2(ThitN)) j . 



Figure I77T1 shows the convergence behavior of L 2 (r(£jv))-error with respect to space 
(in the left subfigure) and time (in the right subfigure) . We observe that the L 2 error 
is of order 0(h 2 ) in space and of 0{At 2 ) in time. 

Figure 17721 shows the convergence behavior of the L 2 (7J 1 )-error with respect to 
space (in the left subfigure) and time (in the right subfigure) . From the left subfigure 
it can be seen that the L 2 (H 1 )-eiTor is of order 0(h) in space. The results in the 
right subfigure indicate that on the meshes used in this experiment the L 2 (iJ 1 )-error 
is dominated by the space error. Note that (very) large timesteps, even At = 1, 
can be used (even for small h), which indicates that the method has good stability 
properties. 

To illustrate the convergence behavior of ii^-errors with respect to time, we con- 
sider an experiment on the shrinking sphere, where the solution is given by u(x, t) = e , 
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(a) (b) 

Fig. 7.1. Convergence w.r.t. L (r(tjv)) norm in Example 1. 





(a) (b) 

Fig. 7.2. Convergence w.r.t. L 2 (i? 1 ) norm in Example 1. 



i.e. the function has maximal smoothness w.r.t. the spatial variable. A simple com- 
putation yields / = || + w-V r M+(div r w)w + A r u = §f-u = 0. The L 2 (iJ 1 )-crrors 
for this example are shown in Figure 17.31 From the left subfigure, we see that again 
the L (if )-error is of order 0(h) in space, just as in the previous case. In the right 
subfigure we observe that the error converges in time with order O(At). 

Example 2. In this example, we consider a surface diffusion problem as in (|2.1|) 
on a moving manifold. The initial manifold is given (as in [5]) by r(0) = { x G 
M 3 | (x\ — x\) 2 + x\ + x\ = 1 }. The velocity field that transports the surface is 

w(x,t) = (O.lxi cos(i), 0.2^2 sin(t), 0.2x3 cos(i)) . 

The initial concentration uq is chosen as uq(x) = 1 + x\XiX^. 

We set At = 0.1 and compute the problem until T = 8. The mesh size of the 
spatial outer mesh is h = 0.125. An approximate surface S h is constructed in the 
same way as in Example 1. In Figure 17741 we show the (aproximated) manifold and 
the discrete solution Uh for different points in time. 

In this problem the total mass M(t) = f r r t \ «(•, t) ds is conserved and equal 
to M(0) = |r(0)| w 13.6083. We check how well the discrete analogon M h (t) = 
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(a) (b) 

Fig. 7.3. Convergence w.r.t. ^(H 1 ) norm for u(x,t) = e t . 

It (t) u h('y t) ds is conserved. In Figure 17751 for different mesh sizes h and a fixed time 
step At = 0.1, this quantity is illustrated. 

From these results we conclude that the method has a satisfactory discrete mass 
conservation property. Note that the errors in the discretization are caused not only 
by the space-time finite element discretization but also by the geometric errors caused 
by the approximation of S by S h . 

8. Conclusions. In this paper we develop a mathematical framework for a new 
Eulerian finite element method for parabolic equations posed on evolving surfaces. 
The discretization method uses space-time elements. The space-time finite element 
method naturally relics on a space-time weak formulation. Such a formulation is 
introduced and shown to be well-posed. The analysis uses a smooth diffcomorphism 
between the space-time manifold and a reference domain. This theoretical framework 
does not allow to treat surfaces that undergo topological changes such as merging or 
splitting. The numerical method, however, can be applied in such situations. Stabiliy 
of the discrete method is derived only for a special case. Numerical experiments 
demonstrate stable behaviour and optimal convergence results also in other cases. 
Extension of the finite element error analysis to more general problems is a topic of 
current research. In this paper, we consider only the case of piecewise linear (in space 
and time) finite elements. The method, however, is directly applicable with higher 
order finite elements. To benefit from the higher order approximation one needs a 
sufficiently accurate approximation of the continuous space-time manifold. 

Acknowledgement. Funding by the German Science Foundation (DFG) through 
the project RE 1461/4-1 and by the Russian Foundation for Basic Research through 
the projects 12-01-91330, 12-01-00283 and 12-01-33084 is acknowledged. The authors 
thank Jorg Grande for his contribution to the implementation of the finite element 
method. 



REFERENCES 

[1] D. Adalsteinsson and J. A. Sethian, Transport and diffusion of material quantities on 
propagating interfaces via level set methods, J. Comput. Phys.. 185 (2003), pp. 271-288. 

[2] B. Alberta, A. Johnson, J. Lewis, M. Raff, K. Roberts, and P. Walter, Molecular 
Biology of the Cell, Garland Science, New York, fouth ed., 2002. 

21 




(e) (f) 

Fig. 7.4. Example 2: solutions for t = 0, 0.4, 2, 4, 6, 



[3] 
[4] 

[5] 

[6] 
[-} 
[8] 
[9] 

[10] 



A.Y. Chernyshenko and M.A. Olshanskii, N on- degenerate Eulerian finite element method 

for solving PDEs on surfaces, Rus. J. Num. Anal. Math. Model., 28 (2013), pp. 101-124. 
A. Demlow and M.A. Olshanskii, An adaptive surf ace finite element method based on volume 

meshes, SIAM J. Numer. Anal., 50 (2012), pp. 1624-1647. 
G. Dziuk, Finite elements for the Beltrami operator on arbitrary surfaces, in Partial differential 

equations and calculus of variations, S. Hildcbrandt and R. Leis, eds., vol. 1357 of Lecture 

Notes in Mathematics, Springer, 1988, pp. 142-155. 
G. Dziuk and C. Elliott, Finite elements on evolving surfaces, IMA J. Numer. Anal., 27 

(2007), pp. 262-292. 
, An Eulerian approach to transport and diffusion on evolving implicit surfaces, Comput 

Visual Sci, 13 (2010), pp. 17-28. 
C. M. Elliott, B. Stinner, V Styles, and R. Welford, Numerical computation of advection 

and diffusion on evolving diffuse interfaces, IMA J. Numer. Anal., 31 (2011), pp. 786—812. 
A. Ern and J.-L. Guermond, Theory and practice of finite elements, Springer, New York, 

2004. 
J. Grande, Finite element methods for parabolic equations on moving surfaces, Preprint 360, 

IGPM RWTH Aachen University, 2013. 

22 




Fig. 7.5. Example 2: Discrete mass conservation 



[11] J. B. Greer, An improvement of a recent Eulerian method for solving PDEs on general ge- 
ometries, J. Sci. Comput., 29 (2008), pp. 321-352. 

[12] S. GROSS and A. ReuSKEN, Numerical Methods for Two-phase Incompressible Flows, Springer, 
Berlin, 2011. 

[13] A.J. James and J. LOWENGRUB, A surfactant- conserving volume- of -fluid method for interfacial 
flows with insoluble surfactant, J. Comp. Phys., 201 (2004), pp. 685-722. 

[14] M. Olshanskii and A. ReuSKEN, A finite element method for surface PDEs: matrix properties, 
Numer. Math., 114 (2009), pp. 491-520. 

[15] M. Olshanskii, A. Reusken, and J. Grande, A finite element method for elliptic equations 
on surfaces, SIAM J. Numer. Anal., 47 (2009), pp. 3339-3358. 

[16] M. Olshanskii, A. Reusken, and X. Xu, A stabilized finite element method for advection- 
diffusion equations on surfaces, IMA J. Numer. Anal., (to appear). 

[17] J. A. Sethian, Theory, algorithms, and applications of level set methods for propagating in- 
terfaces, Acta Numcrica, 5 (1996), pp. 309-395. 

[18] Jian-Jun Xu AND Hong-Kai Zhao, An Eulerian formulation for solving partial differential 
equations along a moving interface, Journal of Scientific Computing, 19 (2003), pp. 573- 
594. 

[19] E. Zeidler, Nonlinear Functional Analysis and its Applications, 11/ A, Springer, New York, 
1990. 



23 



